Synchronization of Micromechanical Oscillators Using Light 
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Synchronization, the emergence of spontaneous order in coupled systems, is of fundamental importance in 
both physical and biological systems. We demonstrate the synchronization of two dissimilar silicon nitride 
micromechanical oscillators, that are spaced apart by a few hundred nanometers and are coupled through optical 
radiation field. The tunability of the optical coupling between the oscillators enables one to externally control 
the dynamics and switch between coupled and individual oscillation states. These results pave a path towards 
reconfigurable massive synchronized oscillator networks. 



Synchronization processes are part of our daily experi- 
ences as they occur widely in nature, for example in fireflies 
colonies 1 , pacemaker cells in the heart 2 , nervous systems 3 
and circadian cycles 4 . Synchronization is also of great tech- 
nological interest since it provides the basis for timing and 
navigation 5 , signal processing 6 , microwave communication 7 , 
and could enable novel computing 8 and memory concepts 9 . At 
the micro and nanoscale, synchronization mechanisms have 
the potential to be integrated with current nanofabrication ca- 
pabilities and to enable scaling up to network sizes. The abil- 
ity to control and manipulate such networks would enable to 
put in practice nonlinear dynamic theories that explain the be- 
haviour of synchronized networks 10 ' 11 . Recent work on cou- 
pled spin torque 7 ' 12 and nanoscale electromechanical oscilla- 
tors (NEMS) 13 ' 14 exhibit synchronized oscillation states. Ma- 
jor challenges with synchronized oscillators on the nanoscale 
are neighbourhood restriction and non-configurable coupling 
which limit the control, the footprint and possible topologies 
of complex oscillator networks 15-17 . Here, we demonstrate the 
synchronization of two dissimilar silicon nitride (S13N4) self- 
sustaining optomechanical oscillators coupled only through 
the optical radiation field as opposed to coupling through a 
structural contact or electrostatic interaction 18,19 . The tunabil- 
ity of the optical coupling between the oscillators enables one 
to externally control the dynamics and switch between cou- 
pled and individual oscillation states. These results pave a 
path towards realizing synchronized micromechanical oscil- 
lators systems connected through optical links. 

Optomechanical oscillators (OMOs) consist of cavity struc- 
tures that support both tightly confined optical modes and 
long-living (high quality factor) mechanical modes 20 ' 21 . These 
modes can be strongly coupled: the cavity optical field leads to 
optical forces acting on the mechanical structure; mechanical 
displacements due to this force in turn affect the cavity opti- 
cal field. Amplification or cooling of the mechanical modes of 
these cavities can be achieved by feeding these cavities with 
a continuous- wave (CW) laser 22 . The mechanical vibration 



(driven by thermal Brownian motion) induces fluctuations of 
the cavity length, which translates into fluctuations of the op- 
tical resonant frequency; for a fixed-frequency driving laser 
this implies that the optical energy stored in the cavity also 
oscillates. Due to a finite cavity optical lifetime, the optical 
field does not respond instantaneously to the mechanical mo- 
tion but instead oscillates with a slight phase lag; as the force 
that the optical field exerts on the mechanical mode is propor- 
tional to the stored optical energy, it will also experience such 
delay. Consequently the optical force will have one compo- 
nent that is in phase with the mechanical displacement, and 
another component that is 90° out of phase. When the laser is 
blue (red) detuned with respect to the optical mode frequency, 
the optical force component that is in phase with the mechani- 
cal vibration induces an optical spring effect that increases (re- 
duces) the mechanical frequency 23 , thus stiffening (softening) 
the mechanical spring. The out of phase component will de- 
crease (increase) the effective mechanical damping, thus am- 
plifying (cooling) the mechanical oscillation. Above a certain 
threshold laser power this optomechanical amplification over- 
comes the intrinsic mechanical damping; the device evolves 
from an optomechanical resonator to a self-sustaining optome- 
chanical oscillator (OMO) 16 . The laser signal fraction that is 
transmitted, or reflected, from the optomechanical cavity be- 
comes deeply modulated at the mechanical frequency of the 
oscillator 20 ' 24 ' 25 . 

Recently it has been predicted that a pair of OMOs could 
synchronize if they are optically coupled as opposed to me- 
chanically coupled 26 ' 27 . Here we experimentally demonstrate 
the synchronization of two optically coupled OMOs [right (R) 
and left (L)] that are fabricated with slightly different dimen- 
sions (i.e. slightly different mechanical frequencies). The 
optical coupling means the mechanical displacement of one 
OMO will lead to a force on the other OMO through the 
optical field. This force is responsible for the effective me- 
chanical coupling between the two OMOs. As the OMOs are 
pumped by a blue-detuned CW laser into self-sustaining os- 
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FIG. 1: Design of the optically coupled optoemchanical oscillators 
(OMOs). (a) Schematic of the device illustrating the mechanical 
mode profile and the optical whispering gallery mode, (b) Scanning 
electron micrograph (SEM) image of the OMOs with chrome heat- 
ing pads for optical tuning by top illumination. (c,d) The symmetric 
(S) and anti- symmetric (AS) coupled optical supermodes. The defor- 
mation illustrates the mechanical mode that is excited by the optical 
field, (e) The dynamics of the coupled OMOs can be approximated 
by a lumped model for two optically coupled damped-driven nonlin- 
ear harmonic oscillators. 



dilations, the R (L) OMO not only experiences the oscilla- 
tion at its natural frequency but also a modulated optical force 
at the L (R) OMO's mechanical frequency. As the coupling 
between the two oscillators is increased, each OMO is eventu- 
ally forced to oscillate at an intermediate frequency between 
their natural frequencies (Qr and Ql), that is, the onset of 
synchronization 27-29 . We observe both the individual free- 
running and synchronized oscillation dynamics by switching 
on and off the purely optical coupling between two OMOs. 

Each individual OMO, shown in figure la,b consists of two 
suspended vertically stacked S13N4 disks. As illustrated in 
figure la, the optical and mechanical modes of such a cav- 
ity are localized around its free-standing edge. The disks are 



fabricated using standard electron-beam lithography followed 
by dry and wet etching steps (see Methods). The two disks 
are 40 um in diameter and 210 nm in thickness, while the 
air gap between them is 190 nm wide. Such a small gap 
and the relative low refractive index of S13N4 (n w 2.0) in- 
duce a strong optical coupling between the top and bottom 
disks. The resonant frequency of the optical modes of the 
stacked disks depend strongly on their separation 30 ; there- 
fore any mechanical vibration that modulates the vertical gap 
width also modulates the optical resonant frequency; a mea- 
sure for the efficiency of this process is the optomechanical 
coupling, defined as g om = duo/dx where uo is the optical fre- 
quency and x is the mechanical mode amplitude 20 ' 30 ' 31 . Our 
device exhibits a large optomechanical coupling, calculated to 
be <7om/27r = 49 GHz/nm (see SI). The mechanical mode that 
couples most strongly to the optical field is also illustrated by 
the deformation of the disks edge in figures la,c which has a 
natural frequency of ^ m /27r « 50.5 MHz. Note that vari- 
ations in the fabrication process lead to different mechanical 
frequencies; indeed we show below that the two cavities are 
not identical and without the optical coupling they oscillate at 
different mechanical frequencies. 

The two OMOs are separated by a distance of d g = (400 ± 
20) nm, minimizing direct mechanical coupling. This gap re- 
sults in evanescent optical coupling between the OMOs when 
their optical resonant frequencies are close. The optical cou- 
pling leads to two optical supermodes spatially spanning both 
OMOs: a symmetric, lower frequency mode &+(£) (figure lc) 
and an anti- symmetric higher frequency mode 6_ (t) (figure 
Id). Their eigenfrequencies are given by uj± = Cj ± k/2, 
where O = (ujl + wr)/2 and ujl (ujr) is the uncoupled op- 
tical resonant frequency of the L (R) OMO and k is the opti- 
cal coupling rate: a reflection of the distance between the two 
cavities. The mechanical modes of each cavity can be approxi- 
mated by a lumped model consisting of two damped harmonic 
oscillators, which are driven by the optical supermode forces, 

Xj^TjXj^Xi = F$ (x R ,x L )/rn$, for 3, k = L,R (1) 

where Xj,Qj,Tj,m^ represent the mechanical displace- 
ment, mechanical resonant frequency, dissipation rate, and ef- 
fective motional mass of each mechanical degree of freedom. 
The optical force is proportional to the optical energy stored 
in the coupled optical modes, which depend both on ^ and 
x L , i.e. F$(x R ,x L ) ex \b±(x R ,x L )\ 2 . Therefore the opti- 
cal field not only drives but also mechanically couples each 
OMO. The nonlinear nature of this driving and coupling force 
form the basis for the onset of synchronization. In a first or- 
der linear approximation when the two OMOs are evenly cou- 
pled (ul = wr), the effective mechanical coupling force be- 
tween the two oscillators is given by F C oup = —kiXj + UqXj 
where kj and Uq are the position and velocity coupling coeffi- 
cients (See supplementary information [SI] for details). These 
coupling coefficients are determined by both the input optical 
power Pi n and laser-cavity detuning A as fcj ex ,, /2 ?2+ A 2)2 
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FIG. 2: Controlling the two OMO system, (a) Schematic of the experimental setup. The pump and probe light are launched together into the 
cavities and are detected separately by photodiodes (PD). An erbium doped fibre amplifier (EDFA) is used to amplify the transmitted signal to 
increase the signal strength, (b) Transmission spectrum of the coupled cavities. The green and orange coloured optical resonances correspond 
to the pump and probe resonances respectively. NT: normalized transmission, (c) Anti-crossing of the optical mode as the relative temperature 
of the L OMO (Tl) and the R OMO (Tr) is changed through varying the tuning laser power. The tuning laser is focused on to the two OMOs 
respectively to obtain the negative and positive relative temperatures, (d) Transmission spectrum of the maximally coupled state indicated by 
the white horizontal line in (c). The red (blue) part of the curve indicates the anti- symmetric (symmetric) optical supermode, k is the optical 
coupling rate. 



and kg oc ((^2)^+A2) 3 m tne unres0AVe d side band limit (op- 
tical damping rate 7 <C Slj) as in our system. Therefore, by 
varying A and P[ n , hence the effective mechanical coupling 
strength, synchronization of the two OMOs can be captured. 

We experimentally demonstrate that the system can be re- 
configured to exhibit either coupled or single OMO dynam- 
ics by controlling the spatial distribution of the optical field 
between the two oscillators. While the distance between the 
two OMOs is fixed (i.e. fixed k), their optical coupling can 
be turned off (on) through increasing (decreasing) the optical 
frequency mismatch 5 = uor — ujl between them. For large 
optical frequency mismatch among the two OMOs (5 ^> k) 
the supermodes reduce to the uncoupled optical modes of the 
individual OMO, (&+,&_) — » (ol,or). This can be read- 
ily seen from the expression of the optical supermodes am- 
plitudes, which are given by linear combinations of the un- 
coupled modes of the left ac(£) and right or(£) cavities: 
b±(t) = a L (t) - a R {t)m/{5 =F (5 2 - k 2 ) 1 ^ 2 ). We tune S ex- 
perimentally using thermo-optic effect, for which the optical 
frequency dependence on temperature can be approximated as 
Wj(Tj) = — g±Tj for j = L,R, where uj^ is the in- 
trinsic optical frequency and g± is the thermal-optic tuning 
efficiency. The thermo-optic tuning is accomplished by focus- 
ing an out-of-plane laser beam with wavelength 1550 nm on 
either OMO (figure 2a). In order to increase the laser absorp- 
tion, we deposit a 200 nm layer of chrome in the centre of 



both OMOs (figure 2a,b). As heat is dissipated in the chrome 
pads, the cavity temperature increases and red shifts the opti- 
cal resonance of the cavity through thermo-optic effect. A sig- 
nature that the optical frequencies of both OMOs are matched 
is given by the almost symmetric resonance dips observed in 
the optical transmission spectrum (figure 2b, 2d), which also 
indicates maximum optical coupling between the cavity opti- 
cal modes. We show experimentally that the coupling of the 
optical modes can be continuously tuned through changing the 
relative cavity temperature as in figure 2c. At AT = we have 
the maximum optical coupling, whereas for AT = ±25 K, the 
relative frequency difference is large (5 ^> k) and the optical 
mode in eq. (1) does not couple the two OMOs. They follow 
the usual single-cavity optomechanical dynamics 20 . 

We first characterize the individual dynamics of the two 
OMOs by switching their optical coupling off. This is 
achieved through increasing the heating laser power such that 
the temperature difference corresponds to the extremities in 
figure 2c. Each cavity is individually excited with a CW laser 
through a tapered optical fibre. As the laser frequency is tuned 
(from a higher to a lower frequency) into the optical reso- 
nance, the radio-frequency (RF) spectrum of the transmitted 
laser signal is detected by a photodiode (PD) and recorded us- 
ing a RF spectrum analyser (RSA). The results revealing the 
single-cavity optomechanical dynamics are shown in figure 
3a,b. The mechanical modes have natural mechanical frequen- 
cies of (f L jR) = (n Li n R )/27r = (50.283,50.219) MHz, 
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FIG. 3: RF spectra of the OMOs and synchronization (a, b) RF power spectra of cavity L (a) and R (b) as a function of laser frequency when 
the coupling is turned off. The horizontal white lines indicate the onset of self-sustaining oscillation. PSD: power spectral density, (c) When 
the coupling is turned on, at an input power Pi n = (1.8 ± 0.2) uW cavities L and R do not synchronize and oscillate close to their natural 
frequencies (see SI), (d) At P; n = (11 ± 1) uW synchronization occurs after the horizontal solid white line. The synchronized frequency 
appears between the two cavities natural frequencies but only appear after a region of unsynchronized oscillation (between the dashed and 
solid white lines), (e) The system oscillate directly in a synchronized state at input optical power P in = (14 ± 1) uW. (f,g,h) Corresponding 
numerical simulations for the OMO system based on the lumped harmonic oscillator model illustrated in fig. Id. NPSD: normalized power 
spectral density. 



and intrinsic quality factors of (Q™ , Q™ ) = (3.4 ± 
0.3,2.3 ± 0.2) x 10 3 . Due to the increased optomechani- 
cal back-action and intracavity optical power the OMOs have 
their frequencies increased (optical spring effect) and ampli- 
tudes grown as the laser is tuned into the optical resonance. 
Above a specific laser-cavity detuning, indicated by the hor- 
izontal white dashed lines on figure 3a,b the intrinsic me- 
chanical losses are completely suppressed by the optomechan- 
ical amplification. At this point the optomechanical resonator 
starts self-sustaining oscillations and becomes an OMO char- 
acterized by sudden linewidth narrowing and oscillation am- 
plitude growth 15, 16,27 . As the laser frequency sweeps away 
from the optical resonance it eventually reaches the point of 
maximum optical power coupled to the cavity; further sweep- 
ing can only reduce the optical power inside the cavity and 
the oscillation vanishes. It is clear from figure 3a,b that each 
cavity has only one mechanical mode in the frequency range 
of interest. Due to the slight difference in geometry, these fre- 
quencies differ by Af = f L - f R = (70.0 ± 0.5) kHz. 

We show the onset of spontaneous synchronization by 
sweeping the CW pump laser across the optical resonance, 
similarly to the single-cavity measurements above only now 
the optical coupling is switched on for coupled dynamics. Us- 
ing the heating laser, we tune the optical coupling to its maxi- 
mum value, indicated by the dashed- white line (Tr — Tl = 0) 
in figure 2c. The laser frequency sweeping is performed at 



various optical power levels corresponding to different ef- 
fective mechanical coupling strength. The optical power 
ranges from slightly above the estimated oscillation thresh- 
old (i.e weaker mechanical coupling for the L and R OMOs, 
Pfa ' ; « (640, 880) nW, up to several times their threshold 
power (i.e. stronger mechanical coupling). At a relative low 
input power, P[ n = (1.8 ± 0.2) uW, the mechanical peaks 
at fn and fa are simultaneously observed on the RF spec- 
trum shown in figure 3c, below the dashed-white line. When 
the laser frequency is closer to the optical resonant frequency, 
more energy is available and the L OMO starts self-sustaining 
oscillation. Since cavity R has a higher oscillation threshold, 
due to its lower mechanical quality factor, it requires more op- 
tical power and only oscillates at a redder detuning; it can be 
noticed from figure 3c that both OMOs oscillate close to their 
natural frequency. Therefore they exhibit asynchronous oscil- 
lations at this lowest power level. At a higher input optical 
power level of P[ n = (11 ± 1) uW, the first oscillation takes 
place at Aujl/2it w —0.10 GHz, and similarly to the case 
shown in figure 3c, the L OMO oscillates first. However, as 
the laser frequency further moves into the optical resonance, 
there is enough energy for both OMOs to start self-sustaining 
oscillations; the two OMOs spontaneously oscillate in unison 
at an intermediate frequency of fs = fts/2ir = 50.37 MHz 
due to the increased effective mechanical coupling, which is 
a clear sign of synchronization. At this time, the output op- 
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FIG. 4: Pump-probe measurement of the individual OMOs oscillation when coupled. The input pump power is Pin = (11 =b 1) uW as in 
fig. 3d. (a,b) Schematic of the pump-probe measurement principle. While the pump laser (green) is symmetrically shared between the two 
OMOs, the probe laser (blue for probing R and red for L) can measure each cavity selectively. (c,d) The uneven probe intensity distribution 
of the cavities, observed by an infrared CCD camera when the pump laser is off. (e) Normalized transmission (NT) spectrum for the probe 
resonances, which correspond to the orange resonances shown in fig. 2b. The red (blue) dashed line corresponds to the probe wavelength 
region for probing the L (R) OMO, as illustrated in (a,b). (f) The red (blue) curve is the L (R) cavity probe transmission RF spectrum, when 
the pump is in the asynchronous region —0.13 < Acjl/2tt < —0.10 GHz shown in fig. 3d; a strong peak at /r is observed but with very 
different amplitude for two probing conditions. The right inset figures show the same curves in linear scale, emphasising the large difference 
between the blue and red curves, (g) Same curves shown in (f) but with the pump laser in the synchronous region Acjl/27t < —0.13 GHz of 
fig. 3d. Here both cavities have similar amplitude at fs, which can be clearly noticed in the linear scale inset. 



tical RF power is increase by more than 5 dB in compari- 
son with the L OMO oscillating only case showing that the 
two OMOs are phase-locked. At an even higher optical in- 
put power, P in = (14 ± 1) uW, the OMOs do not oscillate 
individually, instead they go directly into synchronized oscil- 
lations above the white-dashed line in figure 3e. We confirm 
that the OMOs are indeed synchronized by performing numer- 
ical simulations corresponding to each of the power levels we 
tested. The simulated spectra in figure 3f,g,h exhibit all the 
essential features observed and show good agreement with the 
measured spectra. It also allows us to confirm under which 
conditions the two OMOs are indeed oscillating (see Methods 
and SI). 

To experimentally verify that both structures are indeed os- 
cillating at the synchronized frequency, we probe the mechan- 
ical oscillation of each cavity individually. This demonstrates 
that the single oscillation peaks observed in figures 3d,e are 
not caused by one OMO resonantly driving the other; it also 
verifies that amplitude death of one of the OMOs does not oc- 
cur, a known phenomenon in coupled nonlinear oscillators 32 . 
The transmitted pump laser signal only provides only global 
information of the coupled OMO system, it does not distin- 
guishes the individual contribution from each OMO to the syn- 



chronized signal. To overcome this we used a weak CW probe 
laser, as shown in the setup in figure 2a to excite an optical 
resonant mode that is not strongly coupled between the two 
OMOs; this scheme illustrated in figure 4a,b allows us to se- 
lectively probe the oscillations of the L or R OMO. Figure 
4c, d show the uneven light intensity distribution that can be 
directly observed by capturing the scattered light with an in- 
frared camera. The asymmetric splitting is also evident on the 
probe transmission shown in figure 4e. While these probe opti- 
cal modes exhibit a low optical quality factor (Q op t ~ 4 x 10 4 ) 
that minimizes probe-induced perturbations to the mechanical 
oscillations, the pump laser power and sweep is identical to 
the one used figure 3d. When the L OMO is probed, and the 
pump detuning range is between the dashed and solid lines 
in figure 3d, the probe RF spectrum shows a strong peak at 
/l, which is shown in the red curve in figure 4f. When the 
R OMO is probed, a peak also appears at this frequency, but 
it is 13 dB weaker as shown in the blue curve in figure 4f; a 
weak peak at /r can also be noticed on the blue curve, indicat- 
ing small amplitude oscillatios of the R OMO. These results 
confirm that the oscillation state is very asynchronous in this 
detuning range with the L OMO oscillating at much larger 
amplitude. When the pump laser detuning is above the hor- 
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izontal solid line in figure 3d there is only a single RF peak 
at the synchronized frequency f$ when probing either OMO 
(figure 4g); moreover, they differ in amplitude by less than 0.5 
dB. This shows that both cavities are indeed oscillating with 
similar strength at the synchronized frequency. 

We have demonstrated the onset of synchronization be- 
tween two optomechanical oscillators coupled only through 
the optical radiation field. The ability to control the cou- 
pling strength are promising for realizing oscillator networks 
in which the oscillators can be addressed individually. Fur- 
thermore, established and future micro-photonics techniques 
such as electro-optic and thermo-optic techniques can now be 
extended to switch, filter and phase shift the coupling of these 
oscillators. Here we demonstrated coupling the near field be- 
tween oscillators which can be switched on and off by thermo- 
optical means. In order to achieve long range coupling of 
mechanical oscillators, optical waveguides and optical fibres 
could be used enabling oscillator networks spread over large 
areas only limited by optical waveguide/fibre losses. Optically 
mediated mechanical coupling will also remove the restric- 
tions of neighbourhood while creating 1D/2D/3D mechanical 
oscillator arrays 33 . Using long range, directional and control- 
lable mechanical coupling, synchronized optomechanical sys- 
tems may enable a new class of devices in sensing, signal pro- 
cessing and on-chip non-linear dynamical systems 17 . 

Methods 

Coupled optical cavities A detailed model of this coupled 
optomechanical system is discussed in the SI; here we describe 
a simplified version that capture our system's essential aspects. 
The optical mode amplitudes are approximated by two cou- 
pled harmonic oscillator equations given by 34 , 

ctj = lAjdj — (7/2)07 + mdk + v'Te 5 ?, 

for j,k = L,R,j ^ k 

where Aj = u — Uj(xj,Tj) is the laser-cavity frequency de- 
tuning, oj is the optical frequency of the CW driving laser, 
ujj (xj , Tj ) is the optical frequency of the each uncoupled opti- 
cal mode (clr^l), which depends both on each cavity tempera- 
ture (Tr,l) and mechanical mode amplitude (xr^l). The opti- 
cal dissipation rate is given by 7 and is assumed to be the same 
for both cavities; the coupling rate k measures the evanescent 
field interaction of the two modes and couples the two opti- 
cal cavities. The CW laser drive y/j^Sj, which depends on 
the cavity coupling to the tapered fibre j e , only drives the clr 
mode (-,y%SL = 0), resembling the experimental configura- 
tion shown in figure 2a where the tapered fibre excites only 
the rightmost cavity. Eqs. (2) can be diagonalized to yield the 
coupled optical modes (supermodes) we described in the main 
text. 

Device Fabrication The two 210 nm thick stoichiometric 
SisN4 films are deposited using low-pressure chemical vapour 



deposition (LPCVD). The 190 nm SIO2 layer is deposited by 
plasma-enhanced chemical vapour deposition (PECVD). The 
underlying substrate is a 4 um SIO2 formed by thermal oxi- 
dation of a silicon wafer. The OMOs are defined by electron 
beam lithography which is then patterned by reactive ion etch- 
ing. The heater pads are subsequently defined by photolithog- 
raphy lift-off process. After defining the circular pads with 
lift-off resist, 200 nm of chrome is deposited on the device us- 
ing electron beam evaporation and the residual chrome is lift- 
off afterwards. In order to release the structure, the device is 
immersed in buffered hydrofluoric acid (6:1) for an isotropic 
etch of the SIO2 in between the disks and the substrate layer. 
The device is then dried with a critical point dryer to avoid 
stiction between the two S13N4 disks. 

Experimental setup The schematic for testing the OMO sys- 
tem is illustrated in figure 2a. Two tunable external cavity 
diode lasers are combined using a 3dB directional coupler to 
an optical fiber that is fed into a vacuum probe station. In- 
side the vacuum chamber, the tapered fiber is positioned close 
to the OMO of interest to allow evanecent coupling using a 
micropositioning system. The output light is then splitted by 
a WDM splitter to a New Focus 1811 (125 MHz bandwidth) 
photodetector. Since the power level we use to test for our 
device is low, an erbium doped fibre preamplifer is used to 
amplify the output signal and improve signal-to-noise ratio in 
the detector. The electronic signal from the detector is split 
and fed to an oscilloscope to observe the time waveform and 
to RSA for the frequency spectrum. To obtain the RF map, 
the laser is configured to sweep from the blue side of the reso- 
nance to the red side in a stepwise fashion by applying an ex- 
ternal voltage to the laser cavity piezo-transducer. At each fre- 
quency step, a snap-shot of the RF spectrum is recorded with 
1 kHz resolution bandwidth and 100 Hz video bandwidth. 
Simulation The optical and mechanical quality factors are ob- 
tained by non-linear least square fitting to the measured optical 
and RF spectrum. The measured mechanical and optical fre- 
quencies, and their respective quality factors are fed as param- 
eters to the lumped model described by eqs. (1) and (2). Other 
parameters such as the effective motional mass ra e ff and the 
optomechanical coupling g om are obtained from finite element 
simulations (FEM). From these measured and calculated pa- 
rameters we estimate the threshold power from the expression 
P± = rn Qff Q m uj^(4Q rn Q 3 (g om ) 2 r] c )- 1 , where rj c = J e /(j) 
is the coupling ideality factor. Due to nonlinear dependence 
on most parameters, such as Q and g om the error propagation 
is large in the estimated threshold power. The error propa- 
gation in the threshold estimation is detailed in the SI, where 
we estimate an overall error of dP^/P^ = 35%. To obtain 
the density plots shown in figure 3 we feed the dynamical 
model described in the main text, and detailed in the SI, with 
the measured and calculated parameters. The simulations are 
performed at various power levels close to the experimentally 
measured power. The optimum power matching was chosen 
to match the laser frequencies at which the bifurcations take 
place in the experimental data. 
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SI. EXPERIMENTAL METHODS 
A. Detailed experimental setup 

We measure the optomechanical transduction of the coupled OMOs using the setup shown in figure S 1 . The green (red) line 
indicates the pump (probe) laser path. The probe is only used when taking the pump-probe measurements. The radio frequency 
(RF) spectral maps shown in the main text and in figure SI are obtained with the probe laser off. Both the pump and the probe 
laser are fibre-coupled, tunable, near-infra (IR) lasers (Tunics Reference and Ando AQ4321D). Their optical power is controlled 
using independent variable optical attenuators. The pump and probe light are individually sent to a polarization controller and 
combined with a 50 : 50 directional fiber coupler. A fraction of the power is monitored by a power meter which indicates the 
equivalent input optical power to the system. To prevent the back scattered light from entering the laser, an optical isolator is 
used before feeding the laser into a vacuum probe station (Lakeshore TTPX) operating at a pressure of 10 ~ 5 mT. The light is 
evanescent coupled to the OMOs through a tapered optical fibre waveguide by using a micro positioning system. 

A small portion of the transmitted light (10%) is also monitored by a power meter. The remaining transmitted light is split 
with a wavelength division multiplexing coupler to separate the pump and the probe laser. Since the pump power used is low, 
especially for sub-threshold measurements, the pump light is optionally amplified with a low noise erbium pre-amplifier (EDFA, 
Amonics AEDFA-PL-30) before coupling to a 125 MHz bandwidth photodiode (New Focus 1181). An additional detector 
(Thorlabs PDB150C-AC) can be switched on when the probe measurement is necessary. Half of the detected signal is sent to an 
oscilloscope and the remaining is coupled to a radio-frequency spectrum analyser (RSA, Agilent E4407B). 

The heating light source is provided by another near-IR laser (JDS SWS 16101), operating at 1550 nm, and amplified by a 
high power EDFA (Keopsys KPS-CUS-BT-C-35-PB-lll-FA-FA) that can provide a maximum power of 2 W. The light is sent 
to the microscope optics which focus the light on to the device. Typically, 50 mW of laser power is needed to achieve the desired 
tuning range, details of tuning aspect can be found in section S2 B. 
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FIG. SI: Detailed experimental setup. See SI text for more details. 



B. Measurements 



The RF spectral maps are obtained by detuning the laser from blue to red into the optical resonance in a stepwise fashion, as 
controlled by a voltage applied to laser's external cavity piezo; the laser used has a tuning coefficient of 1.1 GHz/volt. For each 
voltage step, the RF spectrum is recorded. Therefore, the step size determines the vertical resolution of the RF spectra map (see 
Fig. 3 main text) whereas the resolution bandwidth of the RSA determines the horizontal resolution. Here we used a detuning 
step size of 3 MHz and a resolution bandwidth of 1 kHz (100 Hz video bandwidth). This allows us to obtain a high resolution 
map while keeping the data collection time reasonable 20 minutes). 



C. Single and coupled cavities measurement 



The single cavity data are obtained by coupling the tapered fiber either to the L or R OMO. When one OMO is tested, the 
remaining one is heated by the heating laser with high power (~ 50 mW) to ensure that they are completely decoupled. The 
coupled cavity data are obtained by coupling to the R OMO with the tapered fibre. In this case we use the external heating laser 
to fine tune the coupling so that their split spectrum is symmetric. 



D. Pump probe measurement 



The pump probe measurements provide direct evidence for the synchronization of the two OMOs. The individual probe of 
each cavity, as shown in Fig. 4 main text, relies on the asymmetric coupling of one the higher order optical supermodes. This 
asymmetry arises due to their different optical resonant frequency (See section S2 A) which stems from the slight difference in 
the geometry of the two OMOs. This leads to a different mode splitting for the higher and lower order optical modes. In the 
devices we have tested, the majority of them show similar non-identical mode splitting. 

Due to its lowers optical quality factor (Q) and reduced optomechanical coupling g om , the threshold power for self-sustaining 
oscillations 20 ' 35 of the probe resonance is Pth probe ~ 20 mW, which is roughly 20, 000 times larger than the pump resonance 
threshold optical power Pth pump ~ 1 pW. We used a probe power of P p = (20 ± 2) uW, ensuring a low-noise detected probe 
signal without affecting the cavity oscillation dynamics. 
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E. RF time-domain output power 




We record the real-time trace of the output detector signal with an oscilloscope for both the asynchronous state (between 
dashed and solid line in Fig 3c) and the synchronized state (above solid line in Fig 3c). As shown in the RF power spectrum in 
Fig 3d and figure S2 The oscillation output RF power is increased by more than 5 dB as a the two OMO synchronizes. This is 
expected as both OMOs contribute to the total output optical RF power. 



S2. LUMPED MODEL PARAMETERS 
A. Optical and Mechanical modes 



To obtain the optical and mechanical modes of the optomechanical disk cavity we rely on finite element simulations using 
COMSOL®. From these numerical simulations we derive parameters for the lumped model that describes the optomechanical 
dynamics, such as the effective motional mass m e ff, and the optomechanical coupling rate g om . The optical modes are sought 
by solving the Helmholtz vector wave equation with an ansatz E(r, z, <j)) = E(r, z) exp(zra0). In the table SI we show the 
mode radial electric field profile for the lowest order optical transverse-electric (TE) modes. The mechanical displacement field 
is sought by enforcing complete cylindrical symmetry, u(r, </>, z) = u(r, z), the mode profiles are also shown on table SI. From 
the sought eigenmodes, the optomechanical coupling coefficients for the supported optical modes are calculated using boundary 
perturbation theory 36,37 , 



duo ojq 



J(U-n) (Aei 2 \E ■ t\ 2 + Ae^ 1 \D ■ n| 2 ) 



dA 



'° m ~ dx 2 Je\E\ 2 dV 



(SI) 



where the dimensionless displacement field is defined as U = uj max \u\, the relative permitivity differences are given by 
Aei2 =61—62 and Ae^ 1 = 1/ei — l/e2, the unit vectors i and n indicate the tangential and normal components of the vectors. 
The effetive motional mass is calculated as, 



m eff = J p\U\ 2 dV. (S2) 



B. Top illumination thermal tuning 

The coupling between the cavities is controlled by changing their resonant frequencies through the thermo-optic effect. We 
choose to use 200 nm thick chrome pads as the heating element since they absorb 25% of 1550 nm light at normal incidence, 
taking into account its reflectivity. Chrome is also resistant to buffered oxide etch which follows in the fabrication steps. The 
1550 nm laser is amplified with an EDFA, coupled to the imaging microscope and focused on the chrome pads. The heat 
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TABLE SI: Optical and mechanical modes parameters, (a) Geometry of the optomechanical cavity used to calculate the modes and parameters 
shown in the tables. For the optical modes profiles, it is shown the modulus of the radial electric field \E • r|; g om is calculated using 
Eq. (SI), whereas for the mechanical modes it is shown the displacement amplitude \u\ as colors and the deformation represents the normalized 
displacement. 



absorbed by the chrome pads induces a temperature change AT = R t hP a bs, where R t h = dAT/dP a b s ~ 5.2 x 10 3 K/W is the 
simulated effective thermal resistance of our device. Due to thermo-optic effect, the temperature frequency shift rate is given by 
the perturbation expression, 

= Out = fa(r,z)T rel (r,z)\E\ 2 dV 
9th dAT 2n g J\E\ 2 dV ( j 

where < T re \ (r, z) < 1 is the dimensionless relative temperature distribution of the device, a is the material-dependent thermo- 
optic coefficient, and n g is the optical mode group index. If we define the overlap integral T = J SiN \ E\ 2 / J all E\ 2 , Eq. (S3) is 

approximately given by g t h ~ — (T^f*^)uJoGtsiMr/(2rig). In Fig. S3 we show the simulated relative temperature field T re /(r, z), 
at the edge of the disk T re i = T^f*^ ^0.83. From these results we can estimate the top illumination laser power needed to tune 
the cavity's optical frequency by Aujt, 

P abs = _=^_ « _ _ I — I (S4) 



Aujt 2n g ( Aujt\ 

QthRth ~ RtniT^as^T Uo J 



For our device, tuning of SX « 0.2 nm is sufficient to completely decouple the two cavity modes. Using n g « 1.8, asiN 
3 x 10 -5 K _1 , and T « 0.59, Eq. (S4) gives a tuning efficiency gth/Zn ~ —256 MHz/K, therefore a laser power of P 
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P a 6 S /25% w 24 mW is needed to control the optical coupling between the cavities (see section S3 B).This value is in reasonable 
agreement with the experimental power range. 



Chrome heating pad (T = 




Heat sink (T = T ) 



FIG. S3: Thermal tuning of optical resonances. Simulated temperature (AT = T — To) profile of the optical micro cavity. The bottom 
boundary act as a heat reservoir with constant temperature To = 300 K. In the mirroring edge, where the optical modes are localized, the 
temperature is T « 0.83AT 



S3. COUPLED OPTOMECHANICAL OSCILLATORS 
A. Coupled mode equations 

The optical modes a± and a2 of each optical cavity are coupled through the optical near-field. Due to scattering, there is also 
coupling between the clockwise (cw) and counter-clockwise (ccw) optical modes, therefore we need to consider four optical 

. The coupled equations satisfied by these modes are given by 34 ' 38 , 

\ 
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(S5) 



where uj m are optical resonance angular frequencies, j m is total damping rate, k/2 is the inter-cavity optical coupling rate, and 
Vc = le/ilix + 7e) is me coupling ideality factor, where 7 e is the external loss rate (due to the bus waveguide) and 7^ is the 
intrinsic damping rate 39 . 

The system of Eqs. (S5) can be diagonalized exactly, each eigenvector is governed by an equation of the form 



i (mi±) = [-i (Co + (-1)^/2 ± p/2) - 7/2] ^)(i) ffl ^ l(t) ; for m 



1,2, 



(S6) 



where co> = {uj\ + 002)/% 7 = (71 + 72)/2 and £ = ft^l — where J = (71 — 72)/2 + i(co>2 — The original fields 



a cw^,ccw can ^ e recoverec [ f rom t h e eigenvectors through the relation, 



( al w \ / -«6 ( i _) + (£ + W)6( 2 ,-) \ 

JL + (f + ^)&(2,+) 

2£ + (f - ^)&(2 ,-) 

V + (f - ^)&(2,+) y 



x 2 



(S7) 



where &( m ,±) = (fl^ ± Eq. (S7) will be used to calculate the optical transmission function in the section S3 B below. 



B. Steady-state transmission 



To obtain the low-power steady-state optical transmission spectrum, we assume that the laser driving term in Eq. (S5) is 
oscillating at u, i.e., si(t) = s\e lUJt . Eq. (S6) can be written in a rotating frame C( mj ±)(£) = C( mj ±)(t)e lu)t . The resulting 
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equations will be of the form, 



6 (m ,±)= [iA (m , ±) -7/2] b { ^ ±) (±T nSl i^ c ' , form = 1,2, 



(S8) 



where A 



(m,±) 



(cj + (— l) m £/2 ± /3/2) is the laser-cavity frequency detuning for each of the optical supermodes. The 



steady-state solution to (S8) is given by 

km,±) = (TV 



2£ [iA (m>±) - 7/2] 



, for m = 1, 2. 



(S9) 



The driving laser excites directly only the mode a™ , therefore the steady state optical field transmitted through the bus waveguide 
is given by, 

s° 1 ut (oj l ) = s 1 -^T c a c 1 w (S10) 
where the optical field a\ (ui) is given by Eq. (S7). The normalized field transmission, t(uj) = s° ut (uj)/si is given by, 

£ + \an (— iy k 



t = 1 — 1 



2£ 2 ^ 2 V(-l) J /3 + C + 2A + i7 ' (-1)3 -C + 2A + i7 



(SH) 



where A = — co> is the detuning from the average frequency of the two cavities. The normalized power transmission is 
obtained from the relation T(u) = \t(uj)\ 2 . In Fig. S4 we show the transmission T(oS) using the best-fit parameters uj/(2tt) = 
188.442 THz, 7/2tt = 299 MHz, («,/?)/2tt = (1700,298) MHz, and rj c = 0.65. The fit loaded optical quality factor is 
Q = Co/j = 630,000. To obtain the thermal tuned transmission of our device, we use Eq. (SI 1) together with the results 
described in section S2B. The resonant frequency of the cavities, when the top-illumination is on, is given by is given by 
u m (T) = uJmo + 9th AT, where g t h/2n ~ -256 MHz/K (see section S2B). 
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FIG. S4: Optical transmission, (a) Best-fit steady-state normalized optical transmission (red-line), calculated using equation (SI 1), and 
measured transmission spectrum (blue circles). The fit parameters are described in the text, (b) Optical transmission showing the thermal 
tuning of the coupled cavities, the false-color scale indicates the transmission. This map is obtained from (SI 1) using uo\ (T) = cji + g t hAT, 
in good agreement with Fig. 2 in the main text. 



C. Intrinsic mechanical frequencies 



Using the optical read-out of the mechanical motion, as described in section S3 E, we measured the mechanical quality factors 
using a low power optical probe coupled to the lower frequency optical supermodel (right peak in Fig. S4a). The RF spectrum 
(100 averages) showing the two mechanical modes is shown in Fig. S5. The fit parameters are the mechanical frequencies 
and quality factors: (/lJr) = (50.283,50.219) MHz and (Qm L ,Qm R ) = (3.4 ± 0.3,2.3 ± 0.2) x 10 3 . Note that these 
intrinsic frequencies are slightly lower than the OMO self-sustaining oscillation frequency. This is due to the optical spring effect 
explained in the main text. 
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FIG. S5: Mechanical modes RF spectrum, (a) Double-Lorentzian best-fit steady-state normalized RF spectrum (red- line), and measured RF 
spectrum (blue circles). The fit parameters are described in the text. 



D. Mechanical equations and optomechanical coupling 

The mechanical degrees of freedom of each cavity xi, x 2 follows the usual optomechanical equations 20 ' 24 ' 40 ' 41 , 

Si = -riii - nj Xl + -g^- (\a?\ 2 + \a c 2 cw \ 2 ) + if (t), (S12a) 
m K ef J f uo v 7 

x 2 = -T 2 x 2 - n 2 2 x 2 + -g^- (|af | 2 + |af*f ) + if (t), (S12b) 

m eff U 

where f^, 1^, represent the mechanical resonant frequency, dissipation rate, and effective motional mass. 2<t(£) is the ther- 
mal Langevin random force with expectation value (if) = and correlation function (F?(t)F?(t + r)) = 2fcsTm^r^(T), 
where fc^ is the Boltzmann constant and 5(r) is the Dirac delta function. In contrast to the phonon-laser regime 25 , we ignore 
terms which couples, through the mechanical displacement field, the optical modes &(±,i) with &(±,2)5 this is justified because 
k ^> Ql,r- Here we used the optical force as the positive gradient of the energy, this is a convention but must be consistent with 
whether the cavity frequency increases or decreases with increasing mechanical displacement; in our case the optical frequency 
decreases with the mechanical displacement 7 . 

The full optomechanical dynamics is obtained by solving simultaneously Eqs. (S12) and (S5), such dynamics is discussed in 
detail in section S5. It is however instructive to analyze how a prescribed mechanical motion of the two mechanical oscillators is 
read-out through the optical modes (see section S3 E), also how the optical force term in Eqs. (S12) couples to the two of them 
(see section ). 



E. Optical transduction of mechanical oscillations 



To account for the mechanical effect on the optical transmission we first assume that the mechanical motion is independent of 
the optical fields 16 , which is equivalent to ignoring the dynamical back-action. Therefore we can use Eqs. (S6) for the optical 
eigenvectors and simply replace the optical cavity's resonant frequency by uoi — » uoi + g mXi, where X{ is the mechanical 
displacement amplitude for each cavity. The resonant frequency of each eigenmode &( m ,±) will be given by, 

U(i,±)(xi,Xj) =u(xi,Xj) ±£(xi,Xj)/2± P/2, (S13a) 
U(2,±)(xi,Xj) =Q(xi,Xj) ±£(xi,Xj)/2± P/2, (SI 3b) 

where (jj(xi, Xj) = [uJi(xi) +ujj(xj)] /2, £(xi,Xj) = k^/ (1 — [5(xi, Xj)/ n] 2 and 5(xi,Xj) = (ji — 7j)/2 + [ujj(xj) — uji(xi)]. 
Due to the nonlinear £(xi,Xj) dependence on the mechanical displacement amplitudes xi^, The usual analytical approach to 
derive the optomechanical transduction coefficient does not apply 16 . However we can get insight into the problem if we consider 
the strong optical coupling limit, i.e., 5(xi,Xj)/K = g m(xi — xj)/k <C 1 which means that the optical frequency splitting 
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between the cavities is large compared to the mechanically induced frequency shift, therefore £(xi,Xj) ~ k + 0(5 2 /k 2 ). To 
further simplify the analysis we assume that the two cavities share identical optical optical properties, i.e., ui(xi = 0) = 
^2(^2 = 0) = uo and 71 = 72 = 7o- In this case Eq. (SI 3) is approximated by, 

U(m,±)(xi,X 2 ) « Wo(m,±) +W Ol J rX 2 ) (S14) 

where cJo( m ,±) = ^0 + ( — l) m+1 ft/2 ± /3/2. Combined with the above relations, Eq. S6 yields the following equation for the 
optical eigenmodes 6( m> ±), 

&(m,±) = [-i^o(m,±) -i^om^i + ^2) - 7/ 2 ] 6(m,±)(±) m ^ , fori=l,2. (S15) 

The equations above (SI 5) can be formally integrated for a prescribed mechanical motion (xi = AiSm(Qit + 4>i)). The 
homogeneous solutions (si = 0) decay exponentially and does not contribute after the initial transients. To find a particular 
solution satisfying (SI 5) we employ a common approach relying on the Jacobi- Anger expansion 16,27 , 

oo 

exp [ifi, cos(fiit + 00 + i M2 cos(tt 2 t + <f> 2 )] = J2 i P+q J P (Vi)J q (V2)e i{pQl+qQ2)t+i{ * 1+M , (S16) 

p,q— — oo 

where = g om Ai/Qi is the optomechanical modulation depth. Inserting Eq. (S16) in (S15) and solving the resulting equations 
gives, 

2 ^7/2 + i(-Ao (fB , ± )+pni + «n 2 )' (S17) 

where the sum over m, n extends over [— oo, oo], and A ( m5 ±) = uji — ^o(™,±)- From Eq. (S17) we can clearly see the that 
cavity field exhibit tones at combinations of the mechanical frequencies (mQi + nf^) of the two cavities. 

S4. TOY MODEL FOR SYNCHRONIZATION 

In this section we derive an approximate model to describe the essential features of our coupled oscillators. Although we 
develop a first order linear approximation of the two coupled optomechanical oscillators, they constitute an intrinsically a 
nonlinear system, as described in detail elsewhere 27 . 

A. Optically mediated mechanical coupling 

The optical force driving terms in Eqs. (S12) can be written in terms of the diagonal modes &( m ,±) from Eq. (S17) by using 
Eqs. (S7). As in section S3 E, for large optical coupling the terms are only resonant with the driving laser one at a time, therefore 
we can focus on the effect of a particular choice of (ra, ±). To simplify the notation we use A ( m5 ±) = A m and fr( m ,±) = b m 

(i) 

below. We also assume that effective motional mass of the individual oscillators are identical, i.e. = m e ff . The driving force 
in each oscillator is proportional to \b m \ 2 , 

n-P<-™., ^ iP+U p ^i)J q (^)e i ^ + ^ t 2 (slg) 



_ 9om 1 7 1 2 _ QomPjnll Vcl 
^opt — \°m\ — , 

^O(m) 4 ^0(m) 



p,q 



7/2 + i(-A m +g0 2 ) 



which contains both DC terms and oscillatory terms. 

Although our oscillators may exhibit large oscillation amplitude (g m%j > 7), it is instructive to analyze the small amplitude 
dynamics arising for the the driving term in Eq. (SI 8). This treatment is entirely analogous to the one used to derive the 
optomechanical damping and spring effect in uncoupled OMO's 20,41? . For the small amplitude oscillation, the modulation 
parameters are small, i.e., fii = g orn Xi/fti <C 1, therefore the Bessel functions in (SI 8) can be approximated by their small 
argument limit, J n (/i) ~ ^ ( ^) n . We neglect any terms which are quadratic in the /ii ? 2, which also account for summing Eq. 
(SI 8) only over q = 0, ±1 since higher order terms will result in terms which are 0{fi 2 ). The q = terms result in a DC 
component of the force, 

jpU) _ g ° m ( P ^l\r]c\ \ r „ 1Q v 
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The impact of the DC term above is to shift the static equilibirum position of the mechanical oscillators. As a result, the actual 
optical detuning is also shifted, to account for this DC shift we substitute A m — > A f m , where A f m = A m ± g om (xi + #2). 

When p, q = ±1 the resulting terms are quadratic in /ii ? 2 and will be neglected in this first order approximation, therefore 
the lowest order AC terms are given by combinations (p, q) = (0, ±1) and (p, q) = (±1,0). 



■ °Pt AC 



^O(m) 



- £ A, [- cos(^t)/ I (A^, fy) ± sin(%t)/ Q (A^, n,o 



(S20) 



where the functions /j 5 q(A), which correspond to the in-phase (ex sin (Cljt)) and quadrature of phase component (oc cos (Cljt)) 
of the AC force, are given by 



/q(AL,%) 



4(i)A^ 



(( 



2/ 



A£)H 
2A' 



^((2 
((I) 2 



1)2 



AlT 



(( 



i) 2 



AlT 



((I) 4 



A£)+^ 4 ((f) 2 + A£ 



((I) 2 



A /2N 



We can now use the transformations sin(fljt) — > Xj/Aj,cos(fljt) — > xj/(AjQj) and rewrite (S20) as 



rpU) _ 9omPj 
CP 1 AC 



w 0(m) 



■ £ 



-/q(A^,%) 



(S21a) 



(S21b) 



(S22) 



Equation (S22) above shows that for each oscillator the driving force will have a component proportional to its displacement (xj) 
and its velocity (xj). But there are also terms proportional to the displacement and velocity of the opposing OMO; these are the 
terms that couple the two OMOs and form the basis for synchronization in our system. Note that if higher order terms were 
kept in the expansion of Eq. (SI 8), nonlinear terms would appear in Eq. (S22). 

Above we derived the small amplitude form of the optical forces driving our coupled oscillators, we did not use the fact that 
our cavities are in the so-called unresolved sideband regime where the mechanical frequencies are much smaller than the optical 
linewidth, i.e., Qj/j ~ 0.2 <C 1. In this limit, Eqs. (S23) can be written as, 



/q(A^) 

fl( A 'm) 



4(|)A^ 



((f) 2 



+ (Am) 2 ) 

2A' 



3 ' 



Now we can write the (S12) as two coupled harmonic oscillators, 

xx+T^X! ± (^i) 2 xi = 

x 2 + r' 2 x 2 ± (n 2 ) 2 x 2 = 



(A^) 2 

~fax 2 

—k\X\ 



where the modified frequency and damping rate are given by (assuming 5fi| 



k Q x 2 i 
k Q x u 



r' . — r 



Q 



S3) 



(S23a) 
(S23b) 

(S24a) 
(S24b) 

(S25a) 

(S25b) 

(S25c) 
(S25d) 



with/3^ = ^ m P in 7i7/ c i/(m^ ff cj (m)%). 

Therefore, in the small modulation regime (/ii 5 2 <C 1), our system resemble harmonic oscillators in which both the damping 
and frequency are controlled by the optical field; this result is exactly what one would get from uncoupled OMOs. With the 
reduction of the mechanical damping rate for a blue detuned laser. (A^ > 0), these two damped oscillators may undergo a 
bifurcation when the effective damping rate (1^) reverses sign. In this first order approximation there is no additional nonlinearity 
to prevent the oscillations to grow unbound, however it is known that the higher order terms in the force expansion ((S 1 8)) will 
balance the amplitude growth and eventually lead to a stable periodic orbit (limit cycle) 15 ' 16 ' 27 . The optical coupling in our 
system couples the two harmonic oscillators with both amplitude and velocity dependent terms, with coupling strengths fa, kg, 
respectively. The functional dependence of such coupling is the same as the self-induced optical spring and damping rate, as 
given by (S25). 
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B. Approximate Kuramoto model 



In the small amplitude approximation that lead to Eq. (S24), one can also derive slowly-varying phase and amplitude equations 
that describe the dynamics of our system. To accomplish this we assume the following form for our displacement amplitudes, 

x\ (t) = ri (r) exp (r) exp i!7t, (S26a) 
%2 (t) = T2 (t) exp i</>2 (r) exp iQt . (S26b) 

where Q = (Oi + ^)/2 is the average frequency of oscillation. Substituting this ansatz in Eqs. (S24) and assuming that the 
negative damping induced by the optical wave exactly balances the intrinsic viscosity of the oscillators (T' « 0) we obtain the 
following amplitude and phase equations, 

* = + ^ cos ^)) " A ^f^, (S27a) 

* = -"-B + CO < A ^ + ^f^, (S27b) 
A<f> = 2An + - ^ C0S ^ } ^ + 2fcl) + ^ + ^ S ^ A ^ (S27c) 

where A<j) = (pi — fa, R = r\/r2 and Aft = Q[ — ft f 2 . Note that in this approximation the limit cycle has zero amplitude; this 
is the case because we neglected the higher order terms when deriving Eq. (S24). Despite such limitations of this toy model, the 
phase dynamics given by Eq. (S27) enables us to visualize how does the Arnold tongue (| AO| < ki) behaves as we vary the laser 
detuning. We assume that limit cycles of the individual oscillator will have similar amplitude (R = 1) and require A</5 = for a 
synchronized oscillation. When \AQ\ < k\ this conditions can be satisfied and defines an Arnold tongue for this simple model; 
inside the Arnold tongue the system can synchronize (Acj) = 0). In Fig. S6 we show the Arnold tongue plot for our system, in 
Fig. S6a we plot it in the usual way, as function of the coupling coefficient (fcj), whereas in Fig. S6b we plot the tongue as a 
function of the laser detuning by using Eq. (S25). Due to the model simplicity and lack of higher order terms, it does not predict 
the precise values for the synchronization region observed in the experiment ( Fig. S6b), however it does agree qualitatively; for 
higher optical power levels the tongues get wider and allows a given mechanical frequency difference to synchronize at larger 
detuning, as shown in figure 3 of the main text. 



C. Oscillation treshold 



The individual threshold to achieve self-sustaining optomechanical oscillations can be estimated by forcing T' = in Eq. 
(S25)), the resulting threshold condition (A' = 7/2) follows the expression, 



th 



^ m m e ffCJ 4 



(S28) 



where Q is the optical quality factor, Q m is the mechanical quality factor, rj c = 7 e /7 is the ideality coupling factor, it can be 
obtained from the experimental spectrum through its relation to the minimum transmission value 39 . In our experiment rj c = 0.65 
(see Fig. S4). Since the threshold depends nonlinearly on several measured and calculated parameters we propagate the error 
through the usual relation, 




(S29) 



where the 5 before each quantity indicates its standard deviation. The sum of relative errors is dominated by the errors bars of 
g om and Q, using the error bars indicated in the main text we obtain 5P t h/Pth ~ 35%. This is a large deviation but is found to 
impact only the precise optical input power that results in a RF spectral map that matches the experimental data, as shown in the 
main text. 
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Frequency difference (f|-f 2 )/f, % Frequency difference (f i -fi)/f, % 

FIG. S6: Arnold tongue for the simplified Kuramoto model, inside the tongues the system can exhibit synchronized oscillation, (a) Usual 
tongue as a function of the coupling parameter k\. (b) Tongue when k\ is explicitly written in terms of the laser-cavity detuning. The three 
tongues in (b) are obtained with the optical input power values of Pi n = 1, 5, 15 uW, the lower power is the blue whereas the highest power is 
red tongue. / is the average mechanical frequency of the oscillators. 



S5. SYNCHRONIZATION SIMULATION 



A. Simulation approach 

To simulate the synchronization dynamics and obtain the results shown in Fig. S7, we numerically integrate the system of 
equations (S5), including the displacement dependent optical resonant frequencies, i.e. ui^ix) = ^1,2 + 9om x i,2, together 
with the two harmonic oscillator equations (S12). This is accomplished using the NDSolve function in the commercial software 
Mathematical). In the absence of the random thermal noise force in Eq. (S12), it is numerically challenging to capture the 
dynamics before the regenerative oscillation threshold is reached, this is because the steady-state is a static one, i.e., x\p = 0. To 
overcome this issue we add a weak (low-temperature T = 1 K) noise that prevents the dynamics to reach such static equilibrium. 
Since NDSolve is a deterministic solver we include the thermal drive by assigning to F± 2 (t) the outcome of a random variable 
with with expectation value and correlation function given by 

(Ff) = (S30) 

(Ff(t)Ff(t + r)> = 2k B Tmf ff T i 5{r), (S31) 

where fc# is the Boltzmann constant. The discontinuity of this random driving term can lead to instabilities in NDSolve, to 
overcome this we smooth out thenoise term by interpolating the random force with a correlation time t c = (27r/f^)/30. Such 
short correlation time ensures that the noise power spectrum density (PSD) is white within the frequency range of interest.The 
reliability of this approach is confirmed by verifying that for weak pump powers (P <C Pth)> the integrated power spectrum 
density S Xi (Q.) = |xi(^)| 2 satisfy the fluctuation-dissipation theorem 42 . 

^ 2(o) > = ^r^ (o) " n= 2^ (s32) 

A complete analysis of the noise in synchronized systems is beyond the scope of this work, since an accurate numerical noise 
dynamics will require the simulation of the coupled non-linear stochastic dynamics of the optomechanical cavities 43,44 . The 
computational complexity of such systems is also high due to the requirement for slow convergence, first order, fixed time step 
simulation 45-47 . 



20 



B. Simulation results 
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FIG. S7: Numerical simulation of the coupled oscillation dynamics. From a to e: transmission RF spectra, displacement power RF spectra 
of the L and the R OMOs, and the displacement phase diagram of the L and the R OMOs, for input powers at (A)Pi n = 4.9 uW, (B)P* n = 
15.8 uW and (C) P in = 17.9 uW. x L (x R ): displacement of the L and R OMOs. 



The simulation also allows us probe not only the optical transmission PSD, but also the mechanical displacement PSD and 
time series of each OMO. The complete simulation results for the pump laser powers described in the main text are shown in fig. 
S7. The only parameter we adjusted to obtain the maps shown in figures 3 (f,g,h) in the main text and S7 was the optical pump 
power. 

In figures S7A (P in = 4.9 uW), the mechanical power spectrum of the oscillators (fig. S7A(b,c)) shows that for (—0.25 < 
Acj/27t < -0.13 GHz), only L OMO is oscillating; the R OMO is forced to oscillate at the L OMO's frequency but have not yet 
reached its oscillation threshold. This is illustrated by the displacement state space figures shown in fig. S7A(e) for Auj/2tt = 
—0.21 GHz (blue dashed line in fig. S7A(a)), note that \xl | is about 20 times larger than \xr \ . At Auj/2tt = —0.25 GHz, marked 
by the red-dashed line in fig. S7A(a), the situation changes and the R OMO oscillates with larger amplitude (\xr\ w 3.5|#l|) 
but at different frequencies; the result is a Lissajous figure that fills in the whole state space. 

In figures S7B (P in = 15.8 uW), in the asynchronous region, indicated by the blue dashed line, the L OMO oscillates with an 
amplitude roughly 15 times of the R OMO in agreement with the measured RF spectrum and the pump probe measurement. In 
the unified frequency region, for both power levels P in = 15.8 uW and P in — 17.9 uW in fig. S7C, the phase diagram shows 
the two oscillators are synchronized and their amplitude differ less than 20% in agreement to the pump-probe measurements. 
The synchronization phase for figs. S7C(d-e) is roughly = 160°, also all the simulations for our system resulted in phase 
differences close to tt, in agreement with the discussion in 27 that the anti-phase synchronization is a more stable state when the 
oscillations amplitude xl,xr are not identical. 

The criteria for the optimum fitting is the matching of the laser frequency at which the bifurcation occurs, which is sufficient 
for explaining all of the non-linear phenomnena observed. The difference in the simulation power level and experimentally 
measured power level may be due to the variations in etched geometry, film thickness and optical losses. Fitting the entire spectra 
may provide a closer numerical match but it requires a full analysis taking account of of non-linear error propagation, detector 
and spectrum analyzer response function and multidimensional fitting that is beyond the scope of this paper. 



